Finite temperature crossovers near quantum tricritical points in metals 
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We present a renormalization group treatment of quantum tricriticality in metals. Applying a set of flow 
equations derived within the functional renormalization group framework we evaluate the correlation length in 
the quantum critical region of the phase diagram, extending into finite temperatures above the quantum critical 
or tricritical point. We calculate the finite temperature phase boundaries and analyze the crossover behavior 
when the system is tuned between quantum criticality and quantum tricriticality. 
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I. INTRODUCTION 

Quantum critical behavior in metals is a topic of prime in- 
terest for the theory and experiment in the field of condensed 
matter physics^ - — Despite significant effort made over the last 
couple of years, several intriguing puzzles remain unresolved 
by the quickly developing theory of quantum phase transi- 
tions. These include both, fundamental problems regarding 
the correct low energy action to describe quantum critical sys- 
tems, and the physical mechanisms responsible for unusual 
behavior of specific compounds. 

In the conventional scenario of quantum criticality the crit- 
ical temperature T c for a second order transition is driven to 
zero by a tuning parameter, so that a quantum critical point 
emerges. However, it is also possible that below a tricritical 
temperature the transition becomes first order. Then a partic- 
ular brand of quantum criticality arises when the finite tem- 
perature (T) phase boundary terminates at T — with the tri- 
critical point. Strictly speaking, in real systems this scenario 
is hardly feasible, since it requires fine tuning of two non- 
thermal parameters. However, proximity to quantum tricriti- 
cality has been invoked to explain a number of experimentally 
observed properties of metallic compounds.^— 

Although the relevance of quantum tricriticality was em- 
phasized in a number of recent worksr^— a renormalization 
group (RG) study addressing crossovers between critical and 
tricritical behavior, in particular in the quantum-critical, non- 
Fermi liquid regime, is missing. In this paper we present a 
renormalization group analysis for a model system display- 
ing such crossovers. We build upon an earlier work, 14 where 
renormalization group flows of a quantum 6 -model were al- 
ready studied. The main focus of that work was to show that 
order parameter fluctuations may turn first order transitions 
occurring within the bare model into continuous transitions. 

The paper is structured as follows. In Sec. II we present 
a phenomenological picture of the problem to be analyzed, 
and in Sec. Ill we define the corresponding model action. In 
Sec. IV we describe the functional RG framework and derive 
the RG flow equations. In Sec. V we analyze the linearized 
forms of the flow equations obtaining analytical solutions in 
spatial dimensions d > 2. In Sec. VI we present example nu- 
merical results for the correlation length in the quantum crit- 
ical regime, and the shapes of the phase boundaries. In the 



final Sec. VII we lay out our conclusions. 



II. PHENOMENOLOGICAL PICTURE 

Before plunging into the details of the formalism let us con- 
sider the different scenarios which are envisaged in the fol- 
lowing. In Fig.[T](a) we schematically depict a generic phe- 
nomenological phase diagram of a system exhibiting a Gaus- 
sian quantum critical point (QCP). At T — the system can 
be tuned between an ordered and a disordered state by vary- 
ing a non-thermal control parameter r. At finite T, in the dis- 
ordered phase, a crossover region, schematically represented 
here with a line, separates the Fermi liquid and quantum crit- 
ical regimes. For r < 0, the r e -line separates the phases of 
the system at finite T. The classical critical region where non- 
Gaussian fluctuations occur is bounded by the Ginzburg lines 
and vanishes as T — > 0. 

By varying another system parameter the phase diagram 
can be continuously deformed so that the transition becomes 
first order at sufficiently low T, below a tricritical temperature 
r 1 " (see Fig. Q](c)). For a particular value of this parameter 
the tricritical point is located exactly at T — so that the tran- 
sition is second order for all T > but the scaling properties 
of the system for T — » are governed by the quantum tricriti- 
cal point (QTCP) rather than the quantum critical point. This 
special scenario displaying a different scaling behavior is de- 
picted in Fig. [T](b). The purpose of the paper is to compute 
how the crossovers between the scenarios of Fig. Q](a) through 
(b) to Fig.Q](c) occur. Even in scenario (a) a "hidden" tricrit- 
ical point, which would be present at T tn < 0, i.e. the vicinity 
to a first order transition, can affect the scaling behavior at 
higher temperatures. 



III. BARE ACTION 

The conventional description of quantum criticality, which 
we shall rely upon here, invokes the Hertz action. 15 This de- 
scribes a bosonic mode overdamped by particle-hole excita- 
tions across the Fermi level and is applicable under the as- 
sumption that the electronic degrees of freedom may be inte- 
grated out. The original framework of Hertz 15 and Millis 16 
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FIG. 1: (Color online) Different scenarios for generic phase dia- 
grams. Bold solid (dotted) lines mark second (first) order phase tran- 
sitions, thin solid lines indicate the crossover between quantum criti- 
cal and Fermi liquid behavior, while the thin dotted lines indicate the 
boundaries of the Ginzburg region, (a) Schematic phase diagram of a 
system exhibiting a quantum critical point, at which a line of second 
order finite T phase transitions terminates. A crossover region, repre- 
sented here with a line, separates the Fermi liquid and quantum criti- 
cal regimes in the disordered phase, (b) Schematic phase diagram of 
a system exhibiting a quantum tricritical point, (c) Schematic phase 
diagram of a system exhibiting a tricritical point at a finite temper- 
ature. The transition is first order below the tricritical temperature. 
Critical quantum fluctuations are absent. Varying another control 
coupling one can tune between the depicted situations (see the main 
text). 



was recently extended to account for a number of sys- 
tems and phenomena not covered by the original studies. 
These include field-tuned quantum critical points, 17 metam- 
agnetic transitionsr^i^ phase transitions induced by a non- 
equilibrium drivel as well as dimensional crossovers 2 ! and 
quantum criticality involving multiple time scales. 22 We fo- 
cus on the case of discrete symmetry breaking, which can be 
described by a scalar order parameter field <f>, and an action of 
the form 



2^ J (2n 



IPI 



Z-2 



+ Zp z \<t>_ p + U[<p] . (1) 



Here <p p with p = (p, a> n ) is the momentum representation of 
the order parameter field, where a> n = 2nnT with integer n de- 
notes the (bosonic) Matsubara frequencies. Momentum and 



energy units can be chosen such that the prefactors in front 
of jj^jr and p 2 are equal to unity. The action is regularized 
in the ultraviolet by restricting momenta to |p| < A(>. The 
term jjjjpr> where z > 2 is the dynamical exponent, effectively 
accounts for overdamping of the order-parameter fluctuations 
by fermionic excitations across the Fermi surface. The expres- 
sion Eq. (Q~|i is valid for |p|, and ^pr sufficiently small, which 

is the limit relevant for the physical situation considered^ 

The quantity is a local effective potential which we 
expand to sixth order in </>: 



= I dr 



U(tf>(x,T)) 



where 



1/(0) 



Cl2<p 2 + I 



(2) 



(3) 



This ansatz assumes the absence of a field explicitly breaking 
the inversion symmetry. We require a b > to stabilize the 
system at large \(f>\. The bare effective potential yields the well 
known phase diagram, 23 exhibiting a second-order transition 
for a\ > 0, a first order transition for 04 < and a tricritical 
point at «2 = a\ — 0. Within the present model the mass pa- 
rameter cii plays the role of the non-thermal control parameter 
(r) to tune the transition. By varying the quartic coupling 04 
one can deform the phase diagram so that the transition is sec- 
ond or first order at low T (compare Fig.[T|). In conventional 
quantum criticality, the temperature dependence of the coeffi- 
cients d2, a*, as can be neglected, as it leads only to subleading 
corrections. However, it turns out that in a quantum tricritical 
regime the generic quadratic temperature dependence of these 
coefficients yields a dominant contribution to the temperature 
dependence of the correlation length. 



IV. FLOW EQUATIONS 

In the present study we apply the one-particle irreducible 
variant of the functional RGM-— The derivation of the flow 
equations follows Ref. [3, where the present quantum <p 6 
model was applied to analyze the effect of fluctuations on the 
order of quantum phase transitions and the shapes of the finite 
T phase boundaries. The starting point is the exact functional 
evolution equation 2 ^ 



d_ 
OA 



r A M 



1 



d A R A 



2 Tl T< 2 > [cf>] + R A 



(4) 



describing the flow of the effective action r A [<p], which is 
the generating functional for one-particle irreducible vertex 
functions as the infrared cutoff scale A is reduced. Here 
I< 2) [<p] = S 2 r A [<f>]/d(f> 2 , Tr=T 2 W „ / ^0, finally R A denotes 
the cutoff function added to the inverse propagator to cut off 
modes with momentum below the scale A. We implement the 
Litim cutoff 29 



: z(A 2 -p 2 )0(A 2 -p 2 ) . 



(5) 
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For A = Ao the quantity T A is just the bare action, while in 
the limit A — » it converges to the full effective action, i.e., 
the Gibbs free energy. In most cases it is impossible to solve 
Eq. © in the full functional space. One may however cast it 
onto a suitably chosen set of flowing couplings. The approx- 
imation applied here amounts to assuming that the effective 
potential preserves the form given by Eq. ([3]) with flowing 
couplings, while the inverse propagator retains its initial form 



G- l (p,aj n ) = r^[0 = () ] + fl» 

_ \(D, ' 



'Ipl 



-2 



+ Zp 2 + 2a 2 + fl A (p) , (6) 



8,a 4 = \2v d A d 



4 rf (30fl 6 po + 3a 4 ) 
3 2 Z 2 A 4 



■51 



■d fl 6 

ZA 2 



■ 6pod,a 6 



(9) 



c?,fl6 = 16vrfA d 



8 ld (30a 6 po + 3a 4 ) 3 , , , d 30a 6 po + 3a 4 

~3 h + 15 ^ 6 Z 2 A 4 

(10) 

where v d l = 2 d+l n d/2 r(d/2) and t = log(A/A ) < 0. The 
threshold 24 functions l d (S) are determined from 



° 4 Zj^+Ztf+RH9)+Z^5 



(11) 



where (po is the minimum of U(<p). The present study deals 
with situations where the quantum critical point is Gaussian. 
In the following we shall disregard the flow of Z, which is 
equivalent to neglecting the anomalous dimension of the or- 
der parameter field. Such a parameterization of the effective 
action does not capture the anomalous scaling behavior of the 
propagator in the narrow vicinity of the classical second-order 
phase transition at finite temperature. Analogous truncations 
retaining the flow of Z were applied in Refs. fl4H30l The re- 
sults of Ref. 30 indicate that universal aspects of the shape 
of the phase boundary are not affected by non-Gaussian ther- 
mal fluctuations. In the following we shall also neglect the 
renormalization of the factor Z w . The flow of Z w was com- 
puted in Ref. Hoi and shown to be negligible. The present 
truncation reproduces the essential features of the system ex- 
cept for the narrow vicinity of the second order transition at 
T > 0, where nonetheless the shape of the phase boundary 
is described correctly. Note, however, that the approach can 
be adapted to capture the anomalous dimension of the order 
parameter field 30 and to deal with cases where the fixed point 
associated with the quantum critical point is not Gaussian. 31 
The present truncation is analogous to the zeroth order deriva- 
tive expansionj2i2£ which was applied extensively in the con- 
text of classical critical phenomena. The essential quantum 
ingredient present here is the Landau damping term j^t- On 
top of the derivative expansion we performed a polynomial 
expansion of the effective potential. 

Relying on the truncation described above, the flow equa- 
tion Eq. (H]i can be reduced to a set of three ordinary differ- 
ential equations which determine the flow of the couplings a 2 , 
a 4 , fl6 as functions of the cutoff scale A. In practice, whenever 
the effective potential features non-zero minima at +(po, we 
find it more convenient to write the flow equations in terms 
of the variables po = j<Pq, Q4, «6- The coupling a 2 is then 
obtained from 



a 2 = -^paiflA + 3a 6 po) 



(7) 



As long as there exists a non-zero po, the evolution of the 
effective potential is given by the flow equations derived in 
Ref.Q 



4(S) = -Qjfi® ' 



m = - 



i 



n- 1 85 



lUiS), n>2 



(12) 



(13) 



In Eqs. (El SGI) J = S{p) = j^AU'(p) + 2pU"(p)) and the 
threshold functions are evaluated at the value of 6 correspond- 
ing to po, that is, l d n = l„(§)\p =po . The flow equations are illus- 
trated in terms of Feynman diagrams in Fig. [2] 

For po = the flow equations are simplified as the inter- 
action vertices involving an odd number of legs (see Fig. |2| 
vanish. In terms of the couplings a 2 , a 4 , «6 they yield 



d,a 2 - -24v d 



a 4 



1" 

ztf- d l ' 



8,04 = v d A d 



\ 2 ^L)li- 6 Q^l d 



'ZA 2 ) 



ZA 2 1 



(14) 



(15) 



d,a 6 = v d A a 



a 4 



a 6 



12^- ^60 , 
ZA 2 /\ ZA 2 



V 2 



(16) 



V. LINEARIZED FLOW EQUATIONS 



Essential features of the solutions to the equations provided 
in the previous section (the values of the exponents describ- 
ing the system in the quantum-critical regime in particular) 
can be computed analytically by taking only the dominant 
terms in the flow equations into account. Here we analyze the 
linearized version of the flow equations in spatial dimension 
d > 2. We shall use the following rescaled variables: 



2a 2 a 4 
^ZA 2 ' U ^Z^' 



«6 



Z 3 A 6-2</ 



(17) 



We linearize the flow equations around the zero temperature 
Gaussian fixed point 



<9,p = 2v £/ Z-'A d - 2 3+2 



6fl6P0 
6fi!6P0 + «4 



(8) 



8,5 = -25 - 48v d l d (5)u , 



8,u = (d- 4)m - 60v d l d (5)v , 



(18) 



(19) 
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FIG. 2: Feynman diagrams representing the contributions to the flow 
equations for the couplings parameterizing the effective potential. In- 
teractions involving an odd number of legs vanish if the system is in 
the disordered phase. 



d,v = (2d - 6)v . (20) 

All the resulting flow equations involve contributions from 
rescaling. In addition, renormalization of the couplings 6, u by 
the diagrams involving just one interaction vertex is retained. 
After this step, the terms involved are analogous to those an- 
alyzed by Millis, 16 where v = 0, so that only mass renormal- 
ization from the tadpole diagram was considered. Note that 
in the present case of discrete symmetry-breaking, Eq. ( fT8l 
IT9l l20l are applicable both in the symmetric and symmetry- 
broken phases. 

In the above form we still find the flow equations rather 
hard to solve by hand as the threshold function ¥[{6) in- 
volves integrals (see Eqs. (11 1I12I >). Inspired by the seminal 
work by Millis, 16 we additionally expand the function if for 
t = 2 *^f" <k 1 and f » 1. In the regime f <k 1 the flow 
is dominated by quantum physics, while for T » 1 the quan- 
tum contributions to the flow equations are negligible com- 
pared to the classical ones. Similarly to Ref.llg we shall inte- 
grate the flow equations from A = Ao to A = A* determined 
by f(A*) = 1 using the equations valid strictly speaking for 
T <<c 1 . The result will then serve as the initial condition for 
the solution of the flow from A* to A = 0, where we shall in 
turn use the asymptotic flow equations valid in the classical 
sector T » 1. 

Performing the abovementioned expansion of the threshold 
function if we find for f « 1 

4T u 

8,6 * -26 - 48v d -s , (21) 

T d+z-2 



d,u *(d-4)u- 60v d — ; (22) 

T d+z-2 

while for f » 1 

2T 

8,6 « -26 - 48v d — u , (23) 
d 

2T 

8,u = (d- 4)m - 60v d — v . (24) 
d 

The above equations, together with Eq. (1201 ), form a set of lin- 
ear ordinary differential equations of first order. Eq. d20b is 
solved immediately. Plugging the result into the flow of u and 
integrating, we calculate m(A) for A > A* and A < A*. De- 
manding that a r 4 = lirriA^o (24(A) = lirriA^o Z 2 A w(A) > 
we find a condition for the occurrence of a second order tran- 
sition. If a r 4 (T) = at some T > 0, there is a tricritical point 
at this T. We find that the transition is always second order 
at sufficiently high T. It may turn first order as T approaches 
zero. The condition a' 4 > yields 

u >-v (c + C'T^- 2 ^) , (25) 

where uq = u(A = Ao), vo = v(Ao), and C, C are positive 
constants (depending on Ao, Z w , Z, d and z). The explicit ex- 
pressions for these constants are given in the Appendix. If 
the condition for mo given by Eq. d25l l is fulfilled for T — 0, 
it clearly holds for any T > 0. In such case there is a line 
of critical points terminating at T — with a QCP This is 
the scenario (a) in Fig.Q] For vo = one obviously recovers 
mo > and the standard quantum critical behavior. A quantum 
tricritical point exists if uq = = -Cvo, which corresponds 
to the quantum tricritical scenario (b) in Fig. Q] In case the 
condition Eq. d25l l is violated at T — 0, it is still always ful- 
filled at sufficiently high temperatures, namely for T > T m , 
where 

T ai = [-(uo/vd + C)/C'Y Kd+z - 2) ■ (26) 

In this case the transition becomes first order for T < T m , 
which is scenario (c) in Fig.Q] 

We now discuss the flow of 6. Plugging the solution ob- 
tained for w(A) into Eq. ( |2T1 [23) , we integrate Eq. ( |2TI ) from 
Ao to A*. The solution at A* is the initial condition for 
Eq. d23l . In cases where critical behavior occurs, the correla- 
tion length can be extracted as oc lirriA^o A 2 <5. The result 
has the following structure: 

f 2 = D (T) + D\ T (d+z - 2)/z + voD 2 T (2(d+z) - 4)/z , (27) 

where Do depends on 60, uq and vq, while D\ and D2 do not 
depend on 6{). Notice that D\ = D\(uq - M Q n ), where D\ > 0. 
Explicit expressions are quoted in the Appendix. The gener- 
ically quadratic temperature dependence of the bare model 
parameters 60, uq and vo for low T induces a corresponding 
quadratic temperature dependence of Do, D\, D^. While the 
temperature dependence of D\ and D2 is always negligible, 
the quadratic T-dependence of Dq(T) = Do + D' T 2 domi- 
nates in the tricritical regime (see below) and must therefore 
be kept^il 
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Criticality occurs at T — when 5q = 6" such that Do = 0. 
At finite T the behavior of £~ 2 is governed by several terms. 
The quantum critical term £) 1 r (fl,+z_2 ^ z dominates for suffi- 
ciently low T, provided D\ > 0, and reproduces the result for 
£~ 2 (T) derived already by Millis. 16 Compared to that term, the 
quadratic temperature dependence from Dq(T) and the term 
\\)D2T ( ' 2<d+z ^ 4 ^ z coming from the (p 6 interaction are sublead- 
ing corrections. On the other hand, if m = mJ", the coefficient 
D\ becomes zero and the correlation length exhibits a different 
behavior, 



= D' T 2 + v D 2 T 



(2(d+z)-4)/z 



(28) 



The exponent of the second term is 3 forz = 2, d — 3, and 8/3 
for z = 3, d = 3. The latter exponent was obtained already 
from a self-consistent fluctuation resummation. 11 At low tem- 
peratures the "trivial" quadratic term thus dominates the tem- 
perature dependence of £ in the tricritical regime. The rele- 
vance of the quadratic temperature dependences of the bare 
parameters in the tricritical regime was noticed already by 
Misawa et alji 2 *^ However, the exponent for the temperature 
dependence of the order parameter susceptibility obtained by 
these authors for the quantum tricritical regime turned out 
to be identical to that for conventional quantum criticality, 
that is, 3/2 for z - 2. That result, obtained from a self- 
consistent summation of staggered and homogeneous fluctua- 
tions, is clearly at variance with our result^ 

For quantum critical systems close to quantum tricriticality 
(small D\), the correlation length follows a tricritical temper- 
ature dependence above a crossover temperature 



z/(2+z-d) 



(29) 



For the special case of a very small Di, tricritical behavior 
with an exponent 2(flf+ ; ) ~ 4 [ s observed for temperatures above 



Pi 
vqD 2 



zKd+z-2) 



(30) 



Eq. ( f27l > also yields the shape of the phase boundary. From 
the condition £~ 2 = one finds T c ~ \6q - 5"f, where i// = 
i/r qc = d+ z ,_ 2 for quantum criticality, reproducing the result by 
Millis.— This value crosses over to ip - ip m - j when the 
quartic coupling approaches the tricritical value u^ 1 . Only for 
the special case D' Q — one obtains if/ tn = 



2(d+z)-4 



For small 



D' + a crossover between different exponents occurs. 

As already stated, for uq < «[" a classical tricritical point 
exists at T = T m > 0. In Eq. d27l ) this manifests itself by 
a negative value of D\. From the condition £~ 2 = we find 
the shape of the phase boundary above T = T m to follow 
\6 - Sff\ ~ a\T - T tri | + b\T - T tTi \ 2 . This reproduces a MFT 
result for classical tricritical points. 23 

From the solution for (5(A) we also read off the shape of 
the crossover line separating the Fermi liquid and the quan- 
tum critical regimes. The condition for the occurrence of the 
quantum disordered (Fermi liquid) regime 16 is that the corre- 
lation length becomes of the order of the inverse upper cut- 
off before classical scaling is reached, that is, S(A*) > A 2 . 



From this condition we find the shape of the crossover line as 
#0 _ ~ r 2 /% as can also be deduced by a phenomenological 
reasoning. 12 

The linearized flow equations discussed above are not ap- 
plicable in two dimensions. The reason is easily understood 
from the diagrammatic interpretation provided in Fig. [2] The 
tadpole diagram renormalizing the a 2 coefficient exhibits a 
logarithmic infrared singularity for £ — > oo at T > 0. The 
divergence is cured when one accounts for the renormaliza- 
tion of the quartic interaction coupling via the terms of the 
order a 2 . 



VI. NUMERICAL SOLUTIONS TO THE FLOW 
EQUATIONS 



Here we present results from direct numerical solutions to 
the flow equations delivered in Sec. IV. Information concern- 
ing the physical state of the system is read off from the renor- 
malized values of the quantities a 2 , a^, ae in the limit A — > 0, 
a r n = limA^o «n0V). In all numerical computations we put 
Z = 1, Z w = 1, Ao = 1 and a° 6 = «6(Ao) = 1. We omit 
the quadratic T-dependence of a 2 in this section. We analyze 
solutions to the flow equations upon varying the parameters 
a° = cLi(Kq), a° A = «4(Ao) and T. We refer to Ref. 14 for ex- 
ample plots of the numerically computed phase boundaries in- 
cluding cases involving a first order transition. Here we focus 



on the situation where is close to a a} m , a value where a tri- 
critical point exists at T = for some (&. In the analytical cal- 
culation based on the linearized equations one finds the quan- 
tum tricritical point for the choice a° = a ^ - A " 

For d = z = 3, A = a° = 1 this gives af" * -0.0302358. In 
the numerical solution for the full flow equations a value for 
a^' 111 emerges on fine tuning a2 an d a % such that a r 7 w a' 4 « 0, 
where the value for a° tends to come out a bit (by less than 
1%) smaller than the analytical estimate. 

Once the transition line is established, we can fix a® and in- 
vestigate the quantum critical regime spanning over the quan- 
tum critical or tricritical point. In particular, we extract the 
correlation length £ as function of T. In the region of interest 
the system is typically ordered within the bare model (a9 < 0, 
«2 > 0, but (a") 2 > 4a^a^), but becomes disordered in the 
course of the flow, as the order parameter becomes zero for 
some finite A. At this scale we have to switch to the set of 
equations valid in the symmetric phase (see Sec. IV). 

We first show results for pure tricritical scaling of the in- 
verse correlation length £~ 2 as a function of temperature (see 
Fig. [3]). We set d — z — 3 and we choose the parameter 
a 4 ~ a 4 an d a 2 tuned to the zero temperature quantum tri- 
critical point. This corresponds to the situation depicted in 
Fig. 1 (b), where the QTCP is approached from above. Over 
a wide range of temperature we find that the slope e = 2.66 
gives a good fit to the numerical data, corresponding to the an- 
alytically derived value l/i/r 1 " = 8/3. The leveling off at low 
temperatures could be avoided by further fine-tuning a", while 
for higher temperature the slope decreases as non-universal 
features start to play a role towards the cut-off scale. 



j£ 

n(d+z-2) 2 a 6' 
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FIG. 3: (Color online) The logarithm of the inverse correlation length 
versus log(r) in the quantum tricritical regime in d = z = 3. The 
special choice of parameters is a° A = -0.030270317439 and a° 2 = 
0.00018421858. 



If the value of a 4 is increased slightly above a 4 ' n and cr 2 
again tuned to the zero temperature quantum critical point 
the crossover behavior in the correlation length as derived 
in equation ( |27] i appears. Results of the numerical computa- 
tions demonstrating this crossover behavior of the correlation 
length are displayed in Fig. [4] 
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FIG. 4: (Color online) Correlation length versus temperature in the 
quantum critical regime in d = z = 3. The plot parameters are a 4 = 
-0.03 slightly above the tricritical value and a\ = 0.0001809109296 
tuned to the zero temperature transition point. The vertical dashed 
line gives the analytical crossover temperature as obtained from Eq. 
POI l. The fits (straight lines) are discussed in the text. 



At low temperature the slope e — 1.33 fits well to the quan- 
tum critical scaling l/if/ qc = 4/3. For larger temperatures 
the fit gives e = 2.39, which lies about 10% below the value 
expected l/i/r tn . This can be explained with deviations from 
the power law at higher energy as noted in the discussion of 
Fig. |3] On choosing a value of a® closer to the tricritical one 



the crossover temperature decreases. The (low T) region char- 
acterized by the standard quantum critical scaling shrinks and 
tricritical scaling down to T — sets in, which is accompanied 
by e approaching the value 8/3 (see Fig. 3). 

We now turn to the phase boundary. The result computed 
for the quartic coupling a 4 fine-tuned to its tricritical value is 
exhibited in Fig. [5] The shape of the phase boundary follows 
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FIG. 5: (Color online) Transition line in d = z = 3. The quartic 
coupling is chosen close to the tricritical value a° = -0.03027 « 

0,tri 

a 4 . 

the power law T ~ - a^f. In the double logarithmic plot 
we find the slope e « 0.378, which compares well to the value 
ifi** 1 = 3/8 obtained analytically in Sec. V. When the quar- 
tic coupling is increased, the system shows the Hertz-Millis 
scaling with i/r qc =3/4 for low temperature, and at larger tem- 
perature tricritical scaling with tfr m = 3/8 is still visible for 
the right choice of parameters (compare Ref.[l4T). 

VII. SUMMARY 

We have presented a study of crossover behavior occur- 
ring in the vicinity of metallic quantum critical and tricriti- 
cal points. The analysis is based on renormalization group 
flow equations derived within the functional RG framework 
applied to the Hertz action, retaining a <p b term. It comple- 
ments and extends an earlier work (Ref. Il4l) . providing re- 
sults in the quantum critical regime and delivering analytical 
insights. We focused on crossovers occurring in the temper- 
ature dependence of the correlation length and the shape of 
the phase boundary in the quantum critical regime above the 
quantum critical and tricritical points in the phase diagram. 
The linearized form of the flow equations could be solved an- 
alytically in d > 2, yielding exponents for the power-laws 
obeyed by the correlation length and phase boundary. The 
analytical solutions were confirmed and complemented by a 
numerical evaluation of the full integro-differential flow equa- 
tions. In the quantum tricritical regime the power-law contri- 
bution to the inverse correlation length generated by the <p 6 in- 
teraction provides only a subleading correction to the generic 
quadratic temperature dependence induced by the temperature 
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dependence of the coefficients in the bare action. For quantum 
critical systems close to quantum tricriticality, we computed 
the crossover temperature above which tricritical scaling is 
observed. In the situation where a tricritical point occurs at 
T > 0, we obtained an expression for the tricritical temper- 
ature, and recovered the shift exponent corresponding to the 
classical theory of tricriticality. 
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Appendix A: Explicit expressions for the constants in Section V 



For 6(A) one has 



and 



and 



D = (JqAq + A' |w + ^n,AA~J A^ 2 , (A3) 



Dj = A^( M() + v AA^ 



, / 2nZ,, 
B - A ' 



\ z 



A'AA 6 Q - 2d (2nZ a 



D 2 



2 \ Z 



+ 



+B'A^ 2d 



2d-A 2d+z-4 

B[2nZ l \ * !2nZ h ' 
2 — - A \— 



In this Appendix we give explicit expressions for the con- 
stants appearing in the solutions to the flow equations in Sec. 
V. For the equations for u(A) we have introduced 



C = AAl , 



and 



C = An 



B 



2nZ„ 



(Al) 



(A2) 



where 



A = 120vv 



1 



nZ^ (d + z- 2) 1 



A' = 96 A/120 and 



B 



120v rf 
d(d-2)' 



B' = 



96v d 
d(d - 2) 



(A4) 



(A5) 
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